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Abstract 

We summarize the dominant decay modes of the neutral and charged Higgs 
bosons in the Minimal Supersymmetric extension of the Standard Model. 
While two-body decays are in general dominating, the branching ratios for 
three-body decays of the heavy scalar, pseudoscalar and charged Higgs bosons 
can be large below the thresholds if top quarks, W/Z bosons or heavy scalar 
bosons are involved. Analytical expressions have been derived for the partial 
decay widths and the physical implications of these decay modes are discussed. 
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1. Introduction 



The experimental exploration of the electroweak gauge symmetry breaking is one of the 
most important tasks in particle physics. The observation of one or several fundamental 
scalar particles with couplings growing with the masses of the sources would establish the 
Higgs mechanism as the physical basis of the symmetry breaking. The Higgs mechanism 
[0] may be realized in the frame of the Standard Model (SAi) or one of its possible 
extensions among which supersymmetric theories are truly outstanding candidates 0, [| . 
To explore the physical nature of the scalar particles, high-precision measurements of 
their properties are mandatory. To this end precise calculations of the production cross 
sections and the branching ratios of all important decay channels are required 

Since the Higgs couplings to the Standard Model particles are proportional to their 
masses [modulo enhancement /suppression factors in extended models] the most impor- 
tant decay channels are two-body decays to the heaviest particles allowed by phase space. 
Besides these two-body decays, below-threshold three-body decays of the Higgs bosons 
can also be very important. This is well-known in the Standard Model for Higgs decays 
to real and virtual pairs of W/Z bosons in the intermediate mass range ||: although sup- 
pressed by the off-shell propagator and the additional electroweak coupling between the 
W/Z and the fermions, these decay processes are enhanced by the large Higgs couplings 
to the gauge bosons, giving rise to appreciable branching ratios. 

In supersymmetric extensions of the Standard Model, below-threshold decays may 
become even more important, in particular for heavy Higgs decays to virtual and real 
gauge boson pairs, mixed gauge and Higgs boson pairs, as well as top quarks. The analysis 
of three-body decays of the Higgs bosons in the Minimal Supersymmetric Standard Model 
(AtSSAt) has received not much attention in the literature so far. Only recently some 
below-threshold supersymmetric Higgs boson decays have been investigated numerically 
in Ref. 0. We improve on this paper in several aspects: by performing the analysis 
analytically; by completing the ensemble of important decay channels; and last but not 
least, by studying the effects of stop mixing due to non-zero SUSy parameters \i and 
A t . We will perform the analysis in the small and the large limit of tg/5, the ratio of 
the vacuum expectation values, where the results for three-body decays turn out to be 
quite different. Both domains are interesting since they are realized in grand unified 
supersymmetric models with b-r Yukawa coupling unification fTj. 

The paper is organized as follows. In the next Section, we summarize the main features 
of the Higgs sector in the AiSSAi. In Section 3 we discuss the branching ratios in the 
large tg/3 case, with special emphasis on decays for Higgs masses at the edge of the 
SUSy parameter space. In Section 4 we analyze the below-threshold decays of the heavy 
CP-even, CP-odd and charged Higgs bosons for small tg/3 values. The total decay widths 
are summarized in Section 5. Supplementing analytical expressions will be presented in 
the Appendix. 
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2. Physical Set-Up 

In the Minimal Supersymmetric extension of the Standard Model two isodoublets of Higgs 
fields 0] are introduced to provide masses to the up- and down-type fermions. This results 
in a spectrum of a quintet of physical particles: two CP-even neutral scalars h and H, 
one CP-odd neutral (pseudo) scalar A, and a pair of charged scalar particles H ± . Besides 
the four masses, two additional parameters determine the properties of these particles at 
the tree level: the ratio tg(3 of the vacuum expectation values of the two neutral Higgs 
fields and a mixing angle a in the neutral CP-even sector. 

Supersymmetry leads to several relations among the parameters of the Higgs sector, 
and only two of them are in fact independent. If one of the Higgs boson masses [in 
general Ma] and tg(3 are specified, all other masses and the mixing angle a can be derived 
at the tree-level Supersymmetry imposes a strong hierarchical structure on the mass 
spectrum, Mh < Ma 
radiative corrections 



< Mh, My/ < Mh± and Mh < Mz, which however is broken by 



The leading part of the radiative corrections grows as the fourth power of the top 
quark mass mt and the logarithm of the squark mass Ms || . This part is determined by 
the parameter 
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These leading corrections can be summarized in a simple form 
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At the subleading level the radiative corrections are affected by the supersymmetric 
Higgs mass parameter [i and the parameter A t in the soft symmetry breaking interaction^. 
The radiative corrections are positive and they shift the mass of the light neutral Higgs 
boson h upward. The variation of the Higgs boson masses Mh, Mh, Mh± with the 
pseudoscalar mass Ma is shown in Fig. 1 for tg/3 =1.5 and 30 in two mixing scenarios 
[|Kj: (i) (practically^) no mixing fi Ms, A t = 0; and (ii) maximal mixing fi Ms, A t = 
\/QMs- The top quark mass is fixed to m t = 175 GeV and the scalar mass to Ms = 1 



1 We adopt the recent analysis of Ref. JhJ where the full dependence on /i, A t is taken into account and 
where the next-to-leading QCD corrections || are implemented. 

2 A value of fx ~ 100 GeV is not experimentally excluded yet; the partial decay widths we are analyzing, 
change only slightly between /i = and 100 GeV. 
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TeV. If M A is large, the A, H, Higgs bosons are nearly degenerate while the lightest 
Higgs boson h reaches its maximal mass value. Note that for small tg/3 the ordering 
Mu > M H ± > M A holds in the large mass range while M H ± > Mh ~ M A for large tg/3 
instead. The effect of non-zero A t and \x is quite significant since its shifts the maximal 
value of the h mass upward by almost ~ 20 GeV; for large tg/3 values, the relations 
M h ~ M A , Mh — max(M /l ) for M A < max(M h ) and M H ~ M A , M h ~ max(M h ) for 
M A > max(M/J, holding in the case of no mixing, still hold true when the mixing is 
included. 

The couplings of the neutral Higgs bosons to fermions and gauge bosons depend on the 
angles a and /3; they are given in Table 1 with the normalization defined by the SM cou- 
plings 



9H SM f} 



= [V2G F ] 1/2 m f and g H8M vv = 2 [V2G F ] ^ My (4) 



1/2 



The CP-even neutral Higgs bosons h, H share the SM coupling to the massive gauge 
bosons, their couplings to down- (up-) type fermions are enhanced (suppressed) with 
respect to the SM case. If, in the case of large A masses, the mass of the lightest Higgs 
boson is close to its upper limit for a given value of tg/3, the couplings of h to fermions 
and gauge bosons are SM-like while the couplings of the heavy CP-even scalar H are 
suppressed. The pseudoscalar Higgs boson has no tree-level couplings to gauge bosons; 
its couplings to down- (up)-type fermions are (inversely) proportional to tg/3. Radiative 
corrections treated at the level discussed above, are incorporated in the mixing angle a. 
The size of the couplings is shown in Fig. 2a for the mixing scenarios defined before. Here 
again, the mixing has a large impact; however the bulk of the effect consists of shifting 
the h, H masses upward. 





9&UU 


9<S>dd 


9<s>vv 


Hsm 


1 


1 


1 


h 


cos ol j sin/3 


— sin a/ cos /3 


sin(/3 — a) 


H 


sin a j sin f3 


cos a j cos /3 


cos(/3 — a) 


A 


l/tg/3 


tg/3 






Tab. 1: Higgs boson couplings in the MSSM to fermions and gauge bosons relative to the SM. 
Higgs couplings. 

In the limit of large A, H masses, the h mass becomes maximal. In this limit, the 
properties of h approach those of the SM Higgs boson. The heavy neutral scalar Higgs 
boson H decouples from the gauge bosons and the fermionic H couplings approach the 
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corresponding couplings of the pseudoscalar Higgs boson A. For small h, H masses, the 
light Higgs boson h is built-up primarily by the Higgs field H 2 which couples to down- 
type fermions with the strength ~ tg/3, while the heavy Higgs boson H, built-up by Hi, 
couples conversely to up-type fermions. 

For fermions the charged Higgs particles couple to mixtures of scalar and pseudoscalar 
currents, with components proportional to m u ctgf3 and m^tg/3 for the two ± chiralities, 

9H+ud = 

(G F /V2) 1/2 [(1 - j 5 )m u ctg(3 + (1 + j 5 )m d tg(3] (5) 

The couplings of two Higgs bosons with one gauge boson are listed in Table 2. They 
are normalized to the charged/neutral weak couplings 

gw = (V2G F ) 1/2 M W and g z = {V2G F ) 1/2 M Z (6) 

and they come with the sum of the Higgs momenta entering and leaving the vertices. 
Again, radiative corrections are incorporated in the mixing angle a. The magnitude of 
these couplings can be read off Fig. 2a. 
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Tab. 2: The couplings of two Higgs bosons with one gauge boson. They are normalized to the 
weak couplings defined in eq.^) and they come with momenta of the Higgs particles entering 
and leaving the vertices. 

Finally, we summarize the couplings of the three neutral Higgs bosons among them- 
selves which we will need in the subsequent analyses. Unlike the previous cases, the 
radiative corrections are not incorporated in the mixing angle a in total, but additional 
contributions must be taken into account explicitely. Normalized to 

g' z = {V2G F ) l ' 2 M 2 z (7) 

the radiatively corrected self-couplings are to leading order |TT[| 



Afzvih = 2sin2a!sin(/3 + a) — cos2acos(/? + a) + 3— j— — - cos 2 a (8) 

ik/^ sin jj 

^haa = - cos 2(3 cos(/5 + a) + -|— ^ cos 2 (3 

M z sin p 
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These couplings are shown in Fig. 2b for the two mixing scenarios, including subleading 
contributions QT0j| for non-zero [i and A t . 



The value of tg/9 determines to a large extent the decay pattern of the supersymmet- 
ric Higgs bosons. For large values of tg/3 the pattern is simple, a result of the strong 
enhancement of the Higgs couplings to down-type fermions. The neutral Higgs bosons 
will decay into bb and r + r~ pairs, the charged Higgs bosons into tv t pairs below and tb 
pairs above the top-bottom threshold. Only at the edges of the Higgs parameter space 
these simple rules are modified: If h approaches the maximal mass value, the couplings 
become iS.M-like and the decay modes follow the pattern of the Standard Model; if H 
approaches the minimal mass value, it will mainly decay into hh and AA final states. The 
detailed analysis of the large tg/? scenario with special emphasis on 

h -> W*W*,Z*Z* -> 4 fermions (9) 
H -> hh*,AA* -> hbb,Abb (10) 

is presented in Section 3. 

For small values of tg/3 ~ 1 the decay pattern of the heavy neutral Higgs bosons is 



much more complicated pT| , [L2| . The b decays are in general not dominant any more. 
Instead, cascade decays to pairs of light Higgs bosons and mixed pairs of Higgs and 
gauge bosons are important. Moreover, decays to gauge boson pairs play a major role. 
However, for very large masses, the neutral Higgs bosons decay almost exclusively to top 
quark pairs. The decay pattern of the charged Higgs bosons for small tg/3 is similar to 
that at large tg/3 except in the intermediate mass range where cascade decays to light 
Higgs and W bosons are dominant for small tg/3. 

Besides these two-body decays, below-threshold three-body decays of Higgs bosons 
can play an important role. SM. Higgs decays into real and virtual Z pairs provide 
the signature for searching these particles in the intermediate mass range at proton collid- 
ers. The suppression by the off-shell propagator and the additional electroweak coupling 
between the Z boson and the fermions is at least partly compensated by the large Higgs 
coupling to the Z bosons. By the same token, three-body decays of A4SSA4 Higgs 
particles mediated by gauge bosons, heavy Higgs bosons and top quarks, are of physical 
interest, specifically, Fig. 3: 

(i) heavy CP-even Higgs particle H: 

H -> VV* -> Vff [l) (11) 

H -> AZ* -> Aff (12) 

H -> H^W** -> n^ff (13) 

H -> it* -> ibW + (14) 

(ii) CP-odd Higgs boson A: 

A hZ * _> hff (15) 

A -> it* -> ibW + (16) 
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iii) charged Higgs boson H ± : 






H± 


-> hW* 


- hff 




-> AW* 


- Aff 


H± 


-> 6i* -> 


bbW 



(17) 
(18) 
(19) 

The decays H — > VT + £6 can also be mediated by virtual W + H~* and W + W~* interme- 
diate states, and likewise A — > W + ib decays by W + H~* states. However, these contribu- 
tions, given in the Appendix, are very small. [Additional three-body decays, H — > hh* 
or A A* with h*,A* — > 66, could be relevant only for M# a few hundred MeV from its 
minimum value if tg/3 were very large.] Three-body decays of the light CP-even Higgs 
boson h are negligible, with the only exception of h — > W*W* for close to its maximal 
value. 

The decay chains given in the listings above for the heavy Higgs bosons are the dom- 
inant mechanisms in the range of Higgs masses where multi-body decays are relevant. 
This may be exemplified for the first decay chain of A. In a microscopic analysis one 
would consider the chain A — > h*Z* — > (bb)h(ff)z- However, for small values of tg/3, 
the coupling hbb is of order \fGpmb and much smaller than the gauge coupling Zff. 
As a result, the off-shell h* contributions are negligible compared to the off-shell Z* 
contributions below the two-particle threshold. The light Higgs boson h can therefore 
be treated as a stable particle and the decay process can be described as a three-body 
process to a high level of accuracy. This feature enables us to approach the problem of 
below-threshold decays by simple analytical methods. 

In the discussion so far we have disregarded Higgs decays to neutralinos, charginos and 
sfermions. Since we have assumed a common SUSy scale of order Ms ~ 1 TeV, decays 
to sfermions do not play a role in the Higgs mass range of a few hundred GeV which we 
are analyzing. If the chargino/neutralino channels are kinematically open, they affect the 
decay branching ratios in general very strongly ]I2]. Denoting the decay branching ratio 



without neutralino/chargino decays by BR and the branching ratio for neutralino/char- 
gino (and possibly sfermion) decays by BR X , the finally observed branching ratio is given 
by 

BR fm = (1 - BR X ) x BR (20) 

The branching ratios BR X depend on the supersymmetric Higgs parameter /i and the 
gaugino mass parameter M 2 . For large values of fi and M 2 , these decay modes can be 
neglected, yet for parameter sets \i ~ M 2 ~ 0(100) GeV, they could eventually dominate 



over all other decay modes |12| . Since the aim of the present paper is to identify the regions 
in the parameter space where the below-threshold decay modes might be important, we 
will assume, in the subsequent discussion, the decay channels to supersymmetric particles 
to be shut. 
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3. Higgs decay branching ratios: Large tg/3 

For large tg/3 the Higgs couplings to down-type fermions dominate over all other couplings. 
As a result, the decay pattern is in general very simple. The neutral Higgs bosons will 
decay into bb and t + t~ pairs for which the branching ratios are close to ~ 90 % and 
~ 10 %, respectively. In terms of the Fermi decay constant Gp, the partial decay widths 
into fermions are given by [R| 



r($ - //) = n c 



G F M^> 2 2 



(21) 



with p = 3(1) for the CP-even (odd) Higgs bosons; (3 = {1 — Am 2 / Ml) 1 / 2 is the velocity 
of the fermions in the final state, N c the color factor. The couplings <?$// are collected in 
Table 1. For the decay widths into quark pairs, the QCD radiative corrections [|T3J are 
large and must be included. In the limit M$ 3> m g , the 0(a 2 ) corrected decay widths 
read 



T($ -> qq) 



3GpM$ 2 2 



1 + 5.67 



7T 



+ (35.94 - 1.36JV» 



7T 



(22) 



with a s = a s (M$), m 2 = m 2 (M^,); N F = 5 is the number of active quark flavors; all 



0.118 

quark masses evaluated at the scale \i 



226 MeV [which corresponds to the 
The bulk of these QCD corrections can be absorbed into running 



quantities are defined in the MS scheme with re- 
value a s {M\ 



M, 



m ? (M$) = m q (m q ) 



a s {m q 



12/(33-2JV F ) 



1 + cfo 8 (Mg)/7r + c1^(M $ )/tt 2 
1 + cia s (m q )/n + c q 2 a 2 (m q )/n 2 



(23) 



In the case of bottom (charm) quarks, the coefficients are C\ = 1.17(1.01) and c 2 = 
1.50(1.39). For M$ ~ 100 GeV, the b and c quark masses mb(m b ) = 4.23 and m c (m c ) = 
1.23 GeV, as extracted from QCD sum rules Hl5| , have dropped to the effective values^ 
m 6 (M$) ~ 2.9 GeV and m c {M$) ~ 0.62 GeV. 

The charged Higgs particles decay into tv t pairs below and into tb pairs above the 
top-bottom threshold. [Virtual t* — > bW contributions give rise to three-body decays 
for charged Higgs decays — > t*b —>■ Wbb only a few GeV below the two-body decay 
threshold for tb final states.] Neglecting one of the fermion masses in the phase space 
[m u = and m h m t ], the partial decay widths are M 
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-M H ± | V % 



4y/2ir' 



ud 



1 - 



III 
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m 2 d tg 2 (3 + m 2 u ctg 2 (3 - 4 



m u m d 
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(24) 



with Vud being the CKM-like matrix element introduced for quark final states. Also in 
this case the QCD corrections [171 must be included. For light quarks, the bulk of these 



3 A detailed analysis of QCD effects in Higgs decays will be presented in Ref. |lf 
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corrections can be taken into account by using running quark masses; for the top quark, 
the QCD corrections are of order a s /ir and small, since the top quark mass is of the same 
order as the Higgs boson mass. 

The branching ratios for these decays are shown in Fig. 4 for tg/3 = 30. Other de- 
cay modes become important only at the edges of the SUSy parameter space: (i) if h 
approaches the maximum mass for a given value of tg/3, and (ii) if H approaches the 
minimum mass for a given value of tg/3. These two cases are discussed below in detail. 



3.1 The light neutral Higgs boson h 

If the lightest CP-even Higgs boson h approaches its maximum mass for a given value of 
tg/3, the couplings become S.M-like and the decay pattern of the particles is the same as 
for the Higgs particle in the Standard Model. In this part of the parameter space, besides 
decays into charm quarks and gluonsQ, neutral Higgs decays to pairs of real and virtual 
gauge bosons would be observed. Allowing both gauge bosons to be off-shell [important 
only for h masses close to My/ or Mz], the partial decay width reads 

Tlh V*V*\ = 1 ( Ml d Ql M V T V f (Mh - Ql)2 dQlMyTy 

{ ' vr 2 Jo (QI - Ml f + Ml Y\ Jo (QI - Mlf + MyT v 

with Q\ , Q\ being the squared invariant masses of the virtual gauge bosons, My and T v 
their masses and total decay widths. r is given by 



To = ^ls> v sin 2 (/3 - a)X^(Ql Q 2 2 , Ml) 



X(Ql,Q 2 2 ;M 2 h ) + 



12QIQI 

Mi 



(26) 



with 8'y = 2(1) for V = W(Z); A is the two-body phase space function defined by 

\(x, y;x) = (1 — x/z — y/z) 2 — Axy/z 2 (27) 

The branching ratios for the lightest CP-even Higgs boson are shown as a function of 
Mh near its maximal allowed mass value [in practice for Ma = 1 TeV] in Fig. 5a. The 
biggest uncertainty in the prediction of the decay branching ratios is associated with the 
badly known charm quark mass and the QCD coupling constant. To exemplify^ the size 
of the uncertainty we have indicated the error of the branching ratio for h — > cc by the 
shaded band in Fig. 5a. The parameters chosen in the figure are the MS QCD coupling 
a s (M 2 z ) = 0.118 ± 0.006 and m c = 1.23 ± 0.03 GeV for the charm quark mass in the MS 
scheme at the scale of the pole mass as extracted from QCD sum rules ||15|| . The partial 



charm decay width is calculated consistently to 0(a 2 ) with the renormalization scale set 
to the Higgs boson mass [T(J (see also 0]). It is evident from Fig. 5a that the branching 
ratio for h — > cc is very uncertain. Since it is small, the migration of the error to the b 
and r branching ratios is less important, however. 

4 The expression for the gluonic decay width, including the QCD corrections, can be found in Ref. 

5 We refrain from discussing the migration of the uncertainty in the charm sector to the other branching 
ratios which will be displayed in the figure for the central values of the coupling and mass parameters. 
A comprehensive discussion is presented in Ref. p6[ . 
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3.2 The heavy neutral Higgs boson H 

If the mass of the heavy neutral CP-even Higgs boson is close to the minimal value [for 
tg/3 ^> 1, this corresponds to Ma — Mh < 0(100 GeV)], the main decay modes are the 
cascade decays H — > <3><3> with $ = h or A. For real light Higgs bosons, the partial decay 
widths are given by |§ 

^-^otH 1 - 4 !) 1 ' 2 A ™ (28) 

where the radiatively corrected three-boson self-couplings have been given to leading 
order in eqs.(|8|U) and including subleading terms in Fig. 2b. Since Mh — Ma and Xnhh — 
^haa for large values of tg/3, the two branching ratios are equal and close to 50% each. In 
this range of small H masses, the decays H — > $$* with $* — > bb, Fig. 3, can be important 
since the Qbb coupling ~ tg/3 ^> 1 is large. Taking into account only the diagrams where 
the bb final states originate from the decay of h* or A* [the additional channels where h 
or A are emitted from the b lines in H — > bb give negligible contributions since the Hbb 
coupling is very small in this mass range, see Fig. 2a], the Dalitz plot density for this 
three-body decay is given by 

(H -> = Khmj- r=- (29) 



dxidx 2 (1 — X\ — x 2 ) 2 + k$7$ 

where xi j2 = 2Ei >2 /Mh are the scaled energies of the fermions in the final state [which 
we take to be massless], and k$ = Mf/M#,7$ = Y%/M H for $ = h, A. The coefficient 
Kh®® is given by 

js "&G\ M\ 2 2 \2 (nr\\ 

lbn 4 Mh 

The integration of the Dalitz density over the energies x\ and x 2 between the boundaries 

-1 1 K ® 
1 — x 2 — k$ < X\ < 1 

1 - x 2 

< x 2 < 1 - k$ (31) 

can be performed analytically, with the result for Mh ^ 2M$ — T$ sufficiently below the 
thresholds 



Y(H -> <M>*) = IT 



(K$ - 1) ^2 - ^ log K^j + 



1 — 5/«$ / 2k$ — 1 

arctan — , — arctan ■ 



^JAk<$> — 1 V 1/4^$ — 1 y/4K$ — 1 



(32) 



The branching ratios for these two decay modes are shown in Fig. 5b for tg/3 = 30 as a 
function of the pseudoscalar mass and for the two mixing scenarios discussed previously. 
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4. Higgs decay branching ratios: Small tg/3 



In the case tg/3 > 1 the Higgs decay pattern is more complicated because the couplings 
to the light b quarks and r leptons are not enhanced any more and the couplings to the 
heavy top quark and the massive gauge bosons become increasingly important, leading 
in a natural way to a significant admixture of three-body decays. 



4.1 The CP-even Higgs bosons 

Even for small values of tg(3, tg/3 > 1, the light neutral CP-even Higgs boson h decays 
primarily to b quark and r pairs. Only near its maximum mass limit, where the properties 
of h are SM-Mke, other decay modes become non-negligible, as discussed before. The 
branching ratios are shown in Fig. 6a for tg/3 = 1.5 in the cases of "no-mixing" and 
"maximal mixing" introduced above. 

For the heavy neutral CP-even Higgs boson H, besides the usual 2-fermion bb, cc, rr 
and the gg decays, important channels are decays to pairs of light Higgs and gauge bosons 
and mixed pairs of Higgs and gauge bosons. Above the ti threshold, the Higgs boson H 
decays almost exclusively to top quarks. We will concentrate on below-threshold 3-body 
decays involving gauge bosons, Higgs bosons and top quarks, Fig. 3. 

(i) Gauge boson pairs 

In the mass range above the WW and ZZ thresholds, where the HVV couplings are 
not strongly suppressed for small values of tg/3, the partial widths of the heavy Higgs 
particle H for decays into massive gauge bosons may be written 

F(H -> VV) = V F H cos 2 (a - (3){l - Ak v + 12k 2 v ){1 - Ak v ) 1 ' 2 5' v (33) 

327T 

with k v = Ml/Ml and S' v = 2(1) for V = W(Z). Below the threshold for two real 
bosons, the decays to real and virtual VV* pairs are important down to Higgs masses 
M H ~ 130 GeV. The Dalitz plot density for the decay process H -> VV* -> Vff^ is 
given by 

dF [H ^ = Rhvv (1 - xQ(l - x 2 ) + «v(2x t + 2x 2 - 3 + 2k v ) (m) 



dXidx 2 '" ' (1 - X 1 - X 2 ) 2 + Ky~jy 

for all final //' fermions taken massless [i.e. t quark excluded] and summed up. The 
overall normalization factor K HVV is defined as 

3C 2 M 4 

Khw = —^r- cos 2 (/3 - a)M H 5 v (35) 



with 



5 W = 3 and S z = -^-(J^- y s 2 w + ^s 4 w ) (36) 
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The factor 7y = Ty/Mjj accounts for the non-zero width of the gauge boson V. The 
effect of the finite width is important in the threshold region since it guaranties the smooth 
transition from below- to above-threshold decays. Again, Xi are the scaled energies of 
the fermions in the final state. The integration over the energies x\ and x 2 between the 
boundaries in eq. (|3~1"D with k$ replaced by leads to the well-known result of Ref. 
for the partial width [M H < 2M V - T v ] 



t(h -> vv* 



K 



HVV 



1 - 8k v + 20k? 
{Ak v - 1)V2 

1 ~ ^-(2-13^ + 47k 2 



arccos 



6Ky 




6k V + 4K V ) logKy 



(37) 



(ii) Cascade decays 

Potentially interesting decay modes are cascade decays of the heavy Higgs boson H to 
mixed pairs of lighter Higgs bosons h, A and gauge bosons, H — > AZ and H — > if ± iy =F . 
Because of the large value of the top quark mass the radiative corrections will shift the H 
mass upward for a given value of Ma or M#± so that the phase space eventually opens for 
these decay modes. The partial decay widths for the two-body decays above threshold 
are given by the well-known expressions 

T(H^AZ) = sin 2 (/3- a) ^ (M 2 , M 2 ; M 2 H )X(M 2 A , M 2 H ; M 2 ) (38) 

r{H-+H±W*) = ^sin 2 (/3-«)^A^(M| ± ,M^M|)A(M^ ± ,M|;M^)(39) 

with A being the two-body phase space function defined in eq.(^). Below-threshold 
decays of these modes are associated with virtual gauge boson decays to light fermions 
y= t, see later], H — > AZ* — > Aff and H -> if+VF** -> cf. Fig. 3. In analogy to 

the gauge boson pair decays discussed in the previous paragraph, the Dalitz plot densities 
for the two processes may be written 

HF - 3r y2 A/f 4 

-(H^AZ*^Aff) = „ w sm 2 (a - P)M H S z F AZ (x h x 2 ) (40) 



dXidX2 87T 3 

HF 3r^ 2 A/f 4 

(H^H±W**^H±ff') = ^f^ s in 2 (a- (3)M H S w F H±w ( Xl ,x 2 )(Al) 



dxidx 2 87T 3 
where the density function Fij is given by 

F«(*,*,)= ( , t 1 -^ 1 -^-;* (42) 

(1 — X\ — X 2 — Ki + + Kj^fj 

Integrated over the Dalitz plots, the partial widths for decays sufficiently below the thresh- 
olds follow from 



T(H — > AZ* — > Aff) = 6U r w sm 2 (a-(3)M H 5 z G AZ (43) 

an 6 
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T(H -> iZ^W** -> if*//') 



sin 2 (a _ p)M H 5 w G H ± w 



with 



1 

4 



2f-l 



arctan 



(44) 



- A, 



1 — Kj)\/X 



+(X ij - 2Ki) log(/Ci 



IP 



5(1 + «< 



4/Cj Aij 



and 



fmj Top decays 



A 



1 + 2k,; 



(45) 



(46) 



Another class of important decays involves the top quark. Above the top-quark thresh- 
old the partial decay width is given by 

3G F m 2 t 



T(H -> tt) 



9hu M h ft 



(47) 



4\/27r 

with /3f = (1 — 4m 2 /M| r ) 1 / 2 being the velocity of the quarks in the final state; the coupling 
9mt — sin a/ cos/3 has been discussed in the previous section [cf. Table 1 and Fig. 2a]. If 
the mass of the top quark in the Yukawa coupling is evaluated at the scale of the Higgs 
mass mt(Mji), the bulk of the QCD radiative corrections is effectively taken into account. 
Below the tt threshold the final state tbW can be reached through the tt* channel, Fig. 3. 
Since both the W gauge boson mass and the top quark mass must be kept non-zeroQ, the 
Dalitz plot density is more involved than in the previous cases. Using the more convenient 
variables yi t 2 = 1 — £1,2, with x\ t 2 = 2E tt b/MH for the reduced energies of the final top 
and bottom quarks, respectively, and the reduced masses Ki = Mf/M^, one obtains the 
Dalitz plot density 

dr i -, 3G% 9 „,, sin 2 a T , 

<H - W'tb) = -r^mlMl . 2o 2 , (48) 



The amplitude squared T is given by 



647T 3 



sin fiyi + y t Kt 



T = yl(l-yi-y 2 + k w -^K t ) + 2K W (y l y 2 - k w -2K t y l +AK t K W ) 
-K t yiV2 + K t (l- 4« t )(2yi + % + n t ) 



(49) 



When performing the numerical integration over the energies x\,x 2 of the Dalitz plot, 
bounded by 



2(1 - Xi - X 2 + K t + Kb - K>w) + XlX 2 



x\ — An t 



x\ 



4k 6 



< 1 



(50) 



6 Since m^/Af^ <C 1, the b quark mass can be neglected in the amplitude but we shall keep it non-zero 
in the phase space integration where its effect can be important. 
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to obtain the partial decay widths, the non-zero decay width of the virtual particles 
7 4 = rf/Mj^, as well as a non-zero value for kj, = mf/M^, have been included. The 
charge conjugate W + ib final state doubles the partial width. 

Finally, partial widths for H decays into two light Higgs bosons have to be included, 
see eq. (|28|) . For small values of tg/3, one needs to consider only the case where $ = h, A 
are on-shell, since the $66 coupling is very small. 

The branching ratios[] for the main decay channels of the heavy Higgs boson H are 
displayed in Fig. 6b for tg/5 = 1.5. Below the ti threshold, the decay H — > hh is the 
dominant channel, superseded by decays to top quarks above the threshold. This rule 
is only broken for Higgs masses of about 140 GeV where an accidentally small value of 
the Xnhh coupling allows the bb [no mix] and WW* [max mix] decay modes to become 
dominant. Important decay modes in general, below the ti threshold, are decays to pairs 
of gauge bosons, and bb decays. In a restricted range of the H mass, also below-threshold 
AZ* and H^W^* play a non-negligible role. 



4.2 The pseudoscalar Higgs boson A 

Since the pseudoscalar Higgs boson A does not couple to pairs of W and Z bosons, 
the main two-body decay channels are bb and rr decays and, in the mass range between 
about 150 GeV and 350 GeV, cascade decays to hZ final states^. For small tg/5, two-gluon 
decays are competitive in a limited Higgs mass range around 350 GeV. 

(%) Cascade decays 

Above the threshold, the partial width of the cascade decay A — > hZ of the pseu- 
doscalar Higgs boson reads 

T(A -hZ) = ^- cos 2 (/3 - a) ^^(M 2 , Mf; M 2 )A(M 2 , M 2 ; M 2 ) (51) 

with X(x, y; z) being the two-body phase space function defined in eq. (p?D . The branching 
ratio is sizeable for masses between the hZ and the ti thresholds for small values of tg/9. 

For A masses below the threshold, three-body decays A — > hZ* — > hff, mediated by 
virtual Z bosons, can play a role, Fig. 3. In the same notation as before, the Dalitz plot 
density is given by 

dr /If 4 

* .(A - hZ* -> hff) = * w cos 2 (a - (3)M A 5 Z F hz {x x , x 2 ) (52) 



dx\dx2 8n 3 

with all possible Z boson decays to light fermions summed up; the Dalitz plot function 
F is given in eq. (|42|) . Integrating out the energies x\ and x 2 of the fermions, the partial 



7 The Fortran code for all partial decays widths of the SA4 and MSSM Higgs bosons may be obtained 
from djouadi@desy.de, jan.kalinowski@fuw.edu.pl or spira@desy.de. 

8 Since the pseudoscalar Higgs boson is always lighter than H and i? ± in the A4SSA4, the two-body 
decays A — > HZ and A — > iJ ± M /T are kinematically forbidden. 
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width of this decay process reads [Ma < M h + M z — T z ] 



T(A - hZ* -> hff) = * w cos 2 (a - (3)M A 5 z G hz (53) 



with G defined in eg. (|45|) . 

(ii) Top decays 

The three-body process A — > tbW~ , evolving out of the decay 

r(A - «) = (54) 

is predominantly mediated by virtual t quarks, A — > £t* followed by £* — > bW~. Using 
the same notation as for the process H —>■ thW~ we obtain for the Dalitz plot density 

dF -(^^-) ^^,, r ° , (55) 



dxidx 2 647r 3 tg 2 /3 y\ + 

with 

T = yf(l - yi - y2 + k w - Kt) + 2K W (y 1 y 2 - k w ) - n^y^ - 2yi - «t) (56) 

This reduced density approaches the corresponding expression for the scalar CP-even H 
decay H — > tbW~ in the limit Mh,a 3> m 4 as expected from chiral symmetry. The 
non-zero width 7 t = m^/M\ of the top quark and the finite 6-quark mass are taken into 
account when the Dalitz plot density is integrated out to obtain the partial width. 

The numerical results for the two- and three-body decays of the pseudoscalar Higgs 
boson A are summarized in Fig. 6c. The dominant modes are the H — » bb and H —> ti 
decays below the hZ and ti thresholds respectively; in the intermediate mass region, 
Ma = 200 to 300 GeV, the decay H — ► AZ* [which reaches the percent level already at 
Ma = 130 GeV] dominates. Note that this decay mode is more important in the case of 
no mixing. The gluonic decays are significant around the ti threshold. 



4.3 The charged Higgs boson H ± 

Except for the intermediate mass region, the main decay channels of the charged Higgs 
particles are the fermionic tv t and tb modes. The bulk of the radiative QCD corrections 
can be taken into account by using running quark masses. Besides these channels, H + — ► 
cs and cb modes are non-negligible, with the CKM-type suppression of the cb current 
[V c b ~ 0.04] partly compensated by the overall coupling ~ m b of this current. 
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(i) Cascade decays 

If allowed kinematically, charged Higgs bosons can decay into the lightest neutral 
Higgs boson plus a W boson 



hW 



G F cos 2 ((3-a) M, 



w 



8V2- 



TCC 



W 



M H ± 



(Ml M 2 W - M 2 H± )\(M 2 h , M 2 H± - M\ 



2 S 



(57) 



The branching ratio is quite substantial below the top threshold for small values of tg/3 
for which the H ± W T h coupling is not suppressed. The two-body decay of H + to W + A 
is kinematically not allowed. 

There are three below-threshold decays of the charged Higgs boson with substantial 
branching ratios, two of which involve a virtual W boson, H ± — > hW^* — > hff and 
H ± — > AW^* — > Aff, and one involves a virtual t quark, H + — >■ 6t* — > bbW + , Fig. 3. 

The Dalitz plot densities involving the virtual W boson are given by 

dr 



dxida;2 

dr 



dxidx 2 



► hff) 
Aff) 



QG 2 M 4 

* 3 W cos 2 (a-P)M H± F hW (x u x 2 ) (58) 

9G 2 F M^ 

-M H ± F AW (xi 1 x 2 ) 



(59) 



[k, = M 2 / Mfj ± and xi^ being the energies of the // final states], from which one obtains 
the partial decay widths below the thresholds 



T(H + 



> hW + * 
AW + * - 



• hff) 
Aff) 



9G 2 F M^ 



■in 



9G 2 F M^ 



cos 2 (a — (3)M H ± Ghw 
M H ± Gaw 



(60) 
(61) 



The coefficients F and G have been defined in eqs.(f42|, [45l) . 
(it) Top decays 

For the Dalitz plot density of the decay mode H + — > bt* 

3G F X 1/2 (t, b, 



bbW + , evolving out of 



T(H + tb) 



V tb \ z x 



4\/27r M H ± 

"(M 2 ± - m 2 - m 2 ) (m 2 b tg 2 (3 + m 2 t ctg 2 f3) - Am 2 
[with the CKM-type matrix element V t b — 1], we obtain 

dr 



2 2 



(62) 



dxidx 2 
with 



(tf + - 6t* - Wbb) = K H±tb ^-^-^/-w + 2xi + 2x 2 -3 + 2k w 



i-x 2 - n t y + «t7i 



3G 2 m? 1 „, 
32tt 3 tg 2 (3 H 
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(64) 



in the limit where the b quark mass can be neglected; Xi 2 = 2E b i/M H ±. Integrating out 
the energies of the Dalitz plot, the partial width [M#± < m t + — r t ] can be written as 




T{H ± ^ Whh) = l Kn±tb { ^(4K W K t + 3 Kt - 4k w ) log Kw{Kt 1} (65) 

■2 , , , K l ~ 1 5 



+(3/< - 4k* - 3/^ + 1) log 



1 — K W 



4 



-(3« t - KfKw - 2K t K w + Ak w ) + K W 14 — 1 ^ 



The branching ratios of the most important two- and three-body decays, including 
the threshold effects in the numerical analysis, are displayed in Fig. 6d. The branching 
ratio is reduced for the rv search channel of the charged Higgs boson quite significantly; 
indeed, by including of the three-body decays, this decay does not overwhelm all the other 
modes since the H + — > KW* channel as well the channels H — > AW* in the low mass 
range and H + — > bt* in the intermediate mass range, have appreciable branching ratios. 
[Here, as well as for the decay A — > /iZ*, we differ from Ref.Jj].] It is interesting to observe 
that the mixing leads to a qualitative change of the decay pattern in the intermediate 
mass range: in the "no mixing" scenario the hW* decay mode is much more important 
than with mixing for moderately large charged Higgs masses. 



5. Total widths of the Higgs particles 

It is well known that the total widths of the SUSy Higgs particles are in general consid- 
erably smaller than the width of the SM. Higgs particle. 

The large width of the SM. Higgs boson for high masses is due to the decays to 
longitudinal W/Z bosons which grow as GpM^. The absence of the gauge boson couplings 
to A and the decoupling from H for large masses shut this channel. The dominant decay 
modes are built-up by top quarks so that the widths rise only linearly with the Higgs 
masses ~ Gprn? t Mn. Detailed numbers are shown for the two different mixing scenarios 
in Fig. 7. While the width of the SA4 Higgs boson is of order Th sm ~ 200 GeV for a 
Higgs mass M Hsjv[ ~ 750 GeV, the width of the SUSy Higgs particles remain less than 
~ 10 GeV in this mass range for moderate values of tg/3. 

For large tg/3 values, the decay widths of all the five Higgs bosons are determined by 
6-quark finale states and they scale like tg 2 /3 [except for h and H near their maximal 
and minimal mass values, respectively]. The H, A and H + widths therefore become 
experimentally significant, for tg/3 values of order ~ 50 and above and for large Higgs 



masses 21 
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Appendix 

Two of the three-body decays discussed before, H — > tbW and A — > tbW, can be mediated 
by a superposition of different intermediate states. Even though these effects are rather 
small, they should nevertheless be given in this appendix for the sake of completeness. 
Omitting the finite decay widths of the virtual particles (which can easily be implemented) 
and using the notation introduced before, the Dalitz plot density may be written in the 
notation introduced before, 

a) H ^it* + W + H~* + W + W~* -> W + tb : 



dr 



dxidx2 
sin 2 a F 1 



H -> W~tb) 



3G 2 



+ 



64tt 3 
2 sin a sin(/3 



(A.l) 



a) r 



12 



sin 2 /??/ 2 sin/3tg/3 yi(l -yi 
2 sin a cos ((3 — a) T 13 



+ 



sin 2 (/9 — a) T 



22 



-2/2 + %- k+J tg 2 /? (1 - y x 
2 sin(/3 — a) cos(/9 — a) r 23 



V2 + K W - K + 

a) r 



+ 



cos 2 (/3 



33 



sin/3 2/1(1 -yi - y 2 ) tg/3 (1 - y x - V2 + *>w -«+)(l- yi - y 2 ) (1 - y\ - yzj 1 
with the reduced amplitudes squared Ta and the interference terms T^- given by 



Tn = 

Tl2 = 

1^22 — 

Tl3 = 



2/l(l-2/l 
+4K t (2K 



2 
W 



2/2 + - 5/Ct) + 2k w (2/i2/2 - «w0 

- n w n t - k 2 - n w y 1 - 2y 1 K t ) + «t(2j/i - yiy 2 + % + «t) 
(1 - yi - y 2 + k w )(2k w -y\- ym) + 2K t K W (y 2 - yi - 3) 
+«t(2/i + 2/2) (3yi + 2/2 - 1 + 2«t) 
(2/1 + 2/2 - 1 - «w + K t ) \Ak w - + y 2 ) 2 

y 2 (l 2«t) + 3/13/2(1 - 2/i - 2/2 + 3k w - 4«t) 

+2k 2 w {±k w - 2y 2 - 3yi) + 2/t H/ (y 2 + y : - 1) 

+Kt [y 2 (l - y 2 - 2/ct) + yi(l - yi - 2« t ) - 2/^(2/^ + 2/c t + 2yi - 3)] 
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r 



r. 



33 



23 



!/?(% - 1 + Vi + 3y 2 + ««) + ^(-^w - 1 + 2/2 + 3yi + «*) + 2y!y 2 Ki 
—22/13/2 + 2kh/ [j/ 2 («w - 2 + Kf) + ?/i(-k w -2 + n t ) + 2(1 - 

[2/12/2 + % + 2k w -(k w - yi - y 2 )] / «t 
+2/ 2 (l - «t - 3k w - j/i - 3y 2 ) + 2/ 2 (l - «t + «w - S/2 - 3yi) 
+2j/ij/ 2 (l - «t - k w ) ~ *>w + K - l)(yi + 2/2 - 1 + %) 



(A.2) 



b) A—^tt* + W + H~*]\ 



W + tb : 



dxidx 2 



(A -> Wife) 



m 2 M| 
64tt 3 tg 2 /3 



x 



(A.3) 




2/i 2/i(l-2/i-2/2 + «w-«+) (1 - 2/1 - 2/2 + «w - «+) 2 



with the reduced amplitudes squared and the interference term 

Tn = 2/i(l - 2/i - 2/2 + «w - Kt) + 2K W (y 1 y 2 - k w ) - K t {y x y 2 - 2y x - «*) 
Ti2 = (1 - 2/i -1/2 + k w )(2k w - 2/1 - 2/12/2) + «t[(j/i +2/2) 2 - 2/i - 2/2 - 2k w ] 



As expected from chiral symmetry, the expressions approach the corresponding ex- 
pressions for the decay H — > tbW in the limit Mh,a ^ m t- Note that eqs.(A2) and 
(A4) are to be multiplied by a factor of 2 if the charge conjugate processes are taken into 
account. 
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Figure Captions 



Fig. 1: The masses of the MSSM neutral CP-even Higgs bosons h, H and the charged 
Higgs boson as a function of the pseudoscalar A mass for two values of tg/3 =1.5 
(solid lines) and tg/5 = 30 (dashed lines); the long-dashed line presents M H ±. a) 
"no mixing": A t = 0, \x = 100 GeV and Ms = 1 TeV, and b) "maximal mixing": 
A t = y/6M s , n =<^C M s and M s = 1 TeV. [M H ± does not depend on the mixing.] 

Fig. 2a: The coupling parameters of the neutral CP-even Higgs bosons to fermions and 
gauge bosons as functions of the Higgs masses for the two values tg/3 = 1.5 and 
30; no mixing (solid lines) and maximal mixing (dashed lines). The couplings are 
normalized to the SM couplings as defined in Table 1. 

Fig. 2b: The trilinear couplings of the MSSM heavy neutral CP-even Higgs boson, (i) 
Xuhh and (ii) Xhaa as functions of the H mass for the two values tg/3 = 1.5 and 30; 
no mixing (solid lines) and maximal mixing (dashed lines). 

Fig. 3: Main mechanisms for below-threshold three-body MSSM Higgs decays. 

Fig. 4: The branching ratios of the MSSM Higgs bosons h, H, A and H ± as functions 
of their masses for tg/3 = 30; no mixing (solid lines) and maximal mixing (dashed 
lines). For the pseudoscalar and charged Higgs bosons the two mixing scenarios are 
indistinguishable in this case. 

Fig. 5a: The branching ratios of the light CP-even MSSM Higgs boson h as a function of 
its maximum mass value for the two mixing scenarios. The dashed band attached 
to the charm branching ratio, indicates the QCD uncertainties. 

Fig. 5b: The branching ratios of the heavy CP-even MSSM Higgs boson H near its mini- 
mum mass value as a function of Ma for tg/3 = 30 and for the two mixing scenarios. 

Fig. 6a: The branching ratios of the light CP-even MSSM Higgs boson h for tg(3 = 1.5. 

Fig. 6b: The branching ratios of the heavy CP-even MSSM. Higgs boson H for tg/3 = 1.5. 

Fig. 6c: The branching ratios of the CP-odd MSSM Higgs boson A for tg(3 = 1.5; (i) no 
mixing and (ii) maximal mixing. 

Fig. 6d: The branching ratios of the MSSM charged Higgs boson H ± for tg(3 = 1.5; (i) no 
mixing and (ii) maximal mixing. 

Fig. 7: The total widths of the four MSSM Higgs bosons h, H, A and H* 1 as functions 
of the pseudoscalar mass Ma for the two values tg/3 = 1.5 and 30/50 and the two 
scenarios of "no mixing" and "maximal mixing" . 



21 




22 




60 80 100 120 100 200 300 500 700 




Fig. 2a 



23 



- 


! 1 1 


i | 1 | 1 


i 


- 


V \ (i) X(Hhh) 




- 


- 


\ \ 
\ \ 

\ v 
\ \ 


tg/8=30 


- 


— 


1/ A \ 




— 




1/ / \\ 
/ / V 








/i ' X. 






- 


/ i / x ^ _ 

/ / ^ 1 — 




- 




/ 1 / N 

/ N > 
i N 


tg/3=1.5 






no mix 

max mix 


M H [GeV] 






1 1 


, 1,1, 


1 



100 



200 300 



500 700 1000 





1 1 


I I 




(ii) A(HAA) 








— no mix 




Wg|8=30 
\ 

\ 

\ 

V s 


— max mix 








-tgj8=1.5 

y 

✓* 








M H [GeV] 




1 1 


, 1 , 1 , , 



100 



200 



300 
Fig. 2b 



500 700 1000 



24 




25 



.3 
.1 

.03 
.01 

.003 
.001 



— i — i — i — i — i — i — i — — — — — 

— 1 — — 1 — 1 — 1 — — 1 — 1 — 1 — — 1 — h 

= 1 1 bb 1 


i 1 — — 

1 1 1 


I 


1. BR(h) 






- tg,3=30 






Z TT 






— max mix 






— no mix 






gg 




1 1 1 1 


i 1 i i i 1 i i i m 


i4 i! 


l l 



60 80 100 120 140 




1 p 

.3 

.1 b 
.03 

.01 bf 
.003 



.001 



=i=a 



BR(H) 
tg|3=30 



100 
1 
.5 



.01 



100 200 



500 1000 100 
Fig. 4 



tt 



_L 



,M H 
i I i i 



200 300 500700 



n ^ 


1 i i ii M- 








\TV _ 




BR (H*) — 




tg/3=30 : 


cs 


, M H ± 
i i 1 i i i i 



200 



500 1000 



26 



.3 - 



1 — 



.03 



.01 



.003 



i i i i I i i i i I i i i izzr: — i — i — i — r 



BR(h) 
tgj8=1.5 



no mix 



'bb 



TT 




50 



60 



70 



80 



90 



I I I I I I I I I I I I I I I I I I I I -U-UI I I I I -I 



.003 



max mix 




WW/ 
/ M, [GeV] - 

i i i i I i /i i i T i i i i 



50 60 70 80 90 100 110 

Fig. 6a 



28 




100 200 300 500 




29 





30 




100 



200 



300 



500 



BR(A) 
tg/8=1.5 
max mix 



.1 - 



.03 



.01 




500 



31 



"i — i 1 — r 



i — i — r 



=f=i — i — rr 
bt 



BR (H*) 
tg|S=1.5 
no mix 



TV 



X . 
<' \AW 



cs 



cb 




J L 



M H ± 

J i i i 



100 



150 



200 



250 



1 1 1 1 1 1 

TV 


i i i — 1 1 1 — n — = 

1 bt : 


- Aaw* / 


\ BR (H*) - 
\ \ tg/8=1.5 

\ \ max mix 

\ \ 

\ x — 
\ \ — 


V \ / 


\ s — 
\ n — 

\ N — 


A N. 
— / v \ 


\ ^ — 
\ 


/ N \ 
/ \ \ 

/ N \ 
— \ \ 

/ N N 




/ \ 
/ \ 

_ cb / \ 




i i / T^^T^t^ i 


\ \ m h±: 
i \i \\ \ i i i i _ 



100 



150 200 
Fig. 6d 



250 



32 




33 




34 



